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We consider dark matter (DM) that interacts with ordinary matter exclusively through an elec- 
tromagnetic anapole, which is the only allowed electromagnetic form factor for Majorana fermions. 
We show that unlike DM particles with an electric or magnetic dipole moment, anapole dark matter 
particles annihilate exclusively into fermions via purely p-wave interactions, while tree-level annihi- 
lations into photons are forbidden. We calculate the anapole moment needed to produce a thermal 
relic abundance in agreement with cosmological observations, and show that it is consistent with 
current XENON 100 detection limits on the DM- nucleus cross-section for all masses, while lying just 
below the detection threshold for a mass ~ 30 — 40 GeV. 



I. INTRODUCTION 

Roughly 20—25% of the total energy content of the uni- 
verse is in the form of non-baryonic dark matter (DM). 
The exact nature of the dark matter remains a mystery. 
Recently, a variety of authors have explored the possi- 
bility that the dark matter might interact electromag- 
netically with ordinary matter, via an electric or mag- 
netic dipole moment [iT-flij]. Direct detection experi- 
ments strongly constrain such dipole moments for parti- 
cle masses £ 10 GeV. In particular, electric and magnetic 
dipole moments sufficiently small to evade direct detec- 
tion limits cannot provide the correct thermal relic abun- 
dance for the dark matter unless m £ 10 GeV @, 0, EH- 

Although most of the interest in electromagnetic form 
factors has been concentrated on the electric and mag- 
netic dipole moments, Pospelov and ter Veldhuis [l| 
considered several other possible forms of electromag- 
netic coupling to the dark matter, including the electric 
quadrupole moment and the anapole moment. It is the 
latter which we consider in more detail here. The anapole 
moment was first proposed by Zel'dovich jlEj . Unlike the 
electric and magnetic dipole moments, the anapole mo- 
ment has no classical analog, as it does not correspond 
to a multipolar distribution. It is related to the toroidal 
dipole moment, which corresponds to a solenoid with the 
ends joined into a torus, producing an azimuthal mag- 
netic field (for an explanation of the difference between 
the anapole moment and the toroidal dipole moment, see, 
e.g., Ref. 16]). The first experimental measurement of 
an anapole moment in atomic nuclei was noted by Woods 
et al. 

Here, we consider a DM particle that interacts with 
ordinary matter entirely through an electromagnetic 
anapole moment. For Majorana fermions, the anapole is 
the only allowed electromagnetic form factor, making it 
particularly interesting. Anapole dark matter coupling to 
a dark photon was previously considered by Fitzpatrick 
and Zurek |18j], but our model differs from theirs in that 
we have direct coupling to the Standard Model (SM) pho- 
tons, with a correspondingly different Lagrangian. In the 
next section, we begin with the interaction Lagrangian 
for anapole dark matter and explain the general prop- 
erties of the anapole moment. In Sec. IIII1 we use the 



anapole moment operator to derive the total annihilation 
cross-section. We show that it leads, at tree-level, to pure 
p-wave annihilation into fermions. We then use this total 
annihilation cross-section to calculate the relic density of 
the anapole dark matter. In Sec. IIV1 we calculate the 
cross-section for scattering of the anapole dark matter 
particle off of nuclei, and compare this result to the re- 
cent XENON100 limits. We find that, unlike the case 
of the electric or magnetic dipole, anapole dark matter 
with an arbitrary mass can provide the correct thermal 
relic abundance for dark matter and evade current direct 
detection limits. Our results are discussed in Sec. IVj 



II. DARK MATTER ANAPOLE MOMENT 

A Majorana fermion, by definition, is a CPT self- 
conjugate particle. Since the interaction energies for both 
electric and magnetic dipole moments are CPT-odd, a 
Majorana fermion cannot acquire any of these dipole mo- 
ments (although transition electric and magnetic dipole 
moments are still possible 4]). In fact, it has been proven 
that for a massive spin-5 1 Majorana fermion, the only 
electromagnetic form factors allowed are the anapole mo- 
ment and higher multipolcs of it [i~9j . The maximum 
number of these multipoles is 2 5. This theorem has 
been generalized and shown to be valid for any CPT self- 
conjugate massive particle with spin-S 1 , be it fermionic 
or bosonic 20]. For CPT self-conjugate massless par- 
ticles, except the spin-half fermions, no electromagnetic 
couplings are allowed. 

In this article, we consider the dark matter to be a Ma- 
jorana fermion with spin 1/2. According to the theorem 
mentioned above, the only electromagnetic form factor 
allowed is the anapole moment. The interaction operator 
for this anapole moment is of the form xi^ 7 5 X 9 U F^ U 
where \ i s the dark matter and is the electromagnetic 
field strength tensor. (We remark that the correspond- 
ing anapole interaction operator for a Dirac fermion takes 
exactly the same form.) Since this is a dimension-6 op- 
erator, the interaction Lagrangian density is given by 



(1) 



where g is the coupling constant and A is the cut-off 
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TABLE I: Transformation properties of the interaction ener- 
gies for electric dipole moment, magnetic dipole moment and 
anapole moment under C, P, T. 

scale. This interaction operator breaks charge conjuga- 
tion symmetry C and parity symmetry P, but is invariant 
under time-reversal symmetry T. In contrast to the elec- 
tric and magnetic dipole moments which interact with 
external electromagnetic fields, the anapole moment has 
the unique feature that it interacts only with external 
electromagnetic currents J M = d v F iLV . 

The interaction operator for the anapole moment is 
related to but different from that for the toroidal dipole 
moment, which has the form xi 7^ 7 5 Xi d v F^v- One can 
visualize the idea of the toroidal dipole moment as fol- 
lows: a solenoid is folded to join the ends into a torus and 
so the current configuration is such that an azimuthal 
magnetic field is generated. The toroidal dipole vector 
then points in the direction dictated by the curl of the 
magnetic field. Obviously, the toroidal dipole moment 
reduces to the anapole moment in the limit xi = Xi 0- 
In other words, when the incoming and outgoing parti- 
cles are the same, the toroidal dipole moment coincides 
with the anapole moment. This means that we can visu- 
alize the idea of the anapole moment in a similar way as 
that of the toroidal dipole moment. For instance, if neu- 
trinos are Majorana in nature, they can acquire both a 
toroidal dipole moment and an anapole moment formed 
from the various neutrino fields in the mass basis. It has 
been shown that if neutrinos (Dirac or Majorana) indeed 
have a toroidal dipole moment, it will lead to transition 
radiation when the neutrino crosses the interface between 
two media (2l| . 

The transformation property of the anapole moment 
operator under C, P, T is perhaps more transparent in 
the non-relativistic limit, at which the interaction energy 
takes the form 

Hi = -£ff-f, (2) 

where a are the Pauli spin matrices and J — V x B is 
the electromagnetic current. This is consistent with the 
intuitive picture we described above by visualizing the 
folded solenoid, namely we need a non-zero V x B to 
generate an anapole moment. In Table I, we compare 
the transformation properties of the interaction energies 
for the electric dipole moment (a ■ E), magnetic dipole 
moment (a ■ B) and anapole moment (a ■ J) under C, 
P, T in the non-relativistic limit. As one can easily read 
off from Table I, the interaction energy for the anapole 
moment violates C and P individually but preserves T. 
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This is consistent with the transformation properties of 
the anapole interaction operator itself in ((T|) under C, P, 
T. 

A similar but different anapole interaction operator has 
been considered by [l8| and takes the form x 7 M 7 5 X A'. 
If we assume that A'^ is the SM photon, then in or- 
der to maintain gauge invariance, the field x IS re ~ 
quired to transform as x ~ * ^X simultaneously 

when the photon transforms as A 1 ^ — > A'^ — dp£. In- 
deed, the kinetic operator plus the interaction operator, 
X i 7 M dfj_ x + X 7 M 7 5 X A'„ , is invariant under the gauge 
transformations: A' — > A'^ — and x ^ x- So 

the field x acquires a chiral symmetry. However, this is 
impossible unless x 1S massless because the mass term 
m x x X breaks the chiral symmetry. Since a dark matter 
particle must be massive, A 1 ^ cannot be the SM photon. 
Nevertheless, it could still be possible that A 1 ^ represents 
a dark photon which kinetically mixes with the SM pho- 
ton through the operator ef" [H-fH]. 

Therefore, for a fermionic dark matter to couple di- 
rectly to SM photons, the interaction Lagrangian den- 
sity in (JXJ) gives the unique interaction operator for the 
anapole moment. At tree-level, this operator allows for 
the annihilation process XX ~* f f where / is a kine- 
matically allowed SM fermion. (See the Feynman dia- 
gram in Fig. [T]) If we assume that the mass of the 
dark matter m x is smaller than Myy, then the process 
X X W + W~ is not kinematically allowed. The pro- 
cess XX ~~ ^ 77 is kinematically allowed but forbidden. 
(See, for instance, the Feynman diagrams in Fig. [2] and 
Fig. [3 The two similar diagrams with crossed fermion 
lines for Majorana dark matter are not shown.) The rea- 
son for x X ~ 1 7 7 being forbidden at tree-level is that 
d u F /JiV = dfj, {d"A v ) - d 2 A l± , and for on-shell photons, 
both d" ' A v and d 2 A^ are zero. A more intuitive way 
to understand this fact is as follows. The anapole dark 
matter only couples to the external electromagnetic cur- 
rent that generates the electromagnetic fields. But the 
on-shell external photons do not constitute such an elec- 
tromagnetic current. This is in sharp contrast to dark 
matter with electric and magnetic dipole moments, both 
of which allow for the process x X ~ * 7 7 a t tree- level. 
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FIG. 2: t-channel XX ~ > 7 7- 
X * T7 7 




FIG. 3: u-channel XX 
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III. COSMOLOGICAL RELIC ABUNDANCE OF 
ANAPOLE DARK MATTER 

As mentioned in the previous section, for m x < Myp, 
the only annihilation channel for anapole dark matter is 
XX f f where / is a kinematically allowed SM fermion. 
The thermally averaged annihilation cross-section for this 
channel is found to be purely p-wave (see Appendix A for 
the details of the calculation): 
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where a — e 2 /4ir sa 1/137 (note this convention differs 
from that in Ref. [l[) and T is the temperature. 

Using this annihilation cross-section and the method- 
ology of Refs. [Ill [27]], we now derive the thermal relic 
abundance of y. (See also the more recent discussion 
in Refs. [HI, Hlj]). We will confine our attention to 
masses in the range 10 MeV < m x < 80 GeV. The lower 
bound satisfies the requirement that the annihilations of 
the dark matter particles as they become nonrelativis- 
tic in the early universe not heat the photons relative 
to the neutrinos and thereby violate CMB observations 
constraining the number of relativistic degrees of freedom 
; 30j . The upper bound comes from the requirement that 
m x < Myy. Although larger masses are not necessarily 
excluded, we choose to explore that possibility elsewhere. 

We ignore the possibility of an asymmetry between 
particles and antiparticles and also neglect the possibil- 
ity of coannihilations [3l| . Since we have a pure p-wave 
annihilation, we can write the total cross-section as 



( a xx^ff ^1) =<J o 
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where ctq is given by 



cr = 



Ag 2 a m 2 



A 4 



(5) 



and Nj counts the effective number of annihilation chan- 
nels with mass rrif < m x . For each annihilation channel, 
the contribution to Nf is given by the square of the cor- 
responding fermion charge (Q 2 ) multiplied by the color 
factor whenever applicable. For our mass range of inter- 
est, Nf can range from Nf = 1 (for XX ~ * e + e~ only, if 
m x < 100 MeV), up to Nf = 20/3 for m x > m b (anni- 
hilation into 3 charged leptons and 5 quark flavors, with 
each of the latter given by 3 Q 2 for the color factor and 
charge). Following Refs. H,H3| and assuming pure p- 
wave annihilation, we can write the contribution of the 
anapole dark matter to the density in the form 



fl x h 2 = (2.14 x 10 9 ) 



x 2 (GeV)" 1 

1/2 TZ 
g*' Mpi cr 



(6) 



This equation is valid as long as x drops out of thermal 
equilibrium before e + e~ annihilation, which is clearly the 
case for m x > 10 MeV. In Eq. ((6]), O x is the density of % 
relative to the critical density, h is the Hubble parameter 
in units of 100 km sec -1 Mpc" 1 , is the number of rel- 
ativistic degrees of freedom in the universe when x drops 
out of thermal equilibrium, Mpi is the Planck mass, and 
Xf is given by (2(| [27j 
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with g x = 2 being the internal degrees of freedom for 
the Majorana x X pah. If we assume that x X accounts 
for all of the dark matter, then we can substitute the 
observed value of fluMh 2 = 0.11 [32j into Eq. ([6]) and 
use Eqs. © and ([7]) to solve for 00 as a function of m x . 
Then Eq. ([5]) gives the anapole moment, g/A 2 , as a 
function of m x . This is plotted in Fig. 0] Note that the 
step discontinuities are an artifact of the approximation 
used here (in which g* changes sharply as a function of 
temperature and Nf changes sharply as a function of 
mass). Nonetheless, the results shown in Fig. 2] would 
not change significantly with a more detailed calculation. 

In any sensible model, we require g 1, and m x < A. 
This is the case for the range of masses displayed in Fig. 
21 If we set g = 1 for naturalness, then as m x is varied 
from 10 MeV to 80 GeV, A varies from 2.2 GeV to 340 
GeV. 

Finally, while it is not shown in Fig. 21 we remark 



that anapole dark matter with mass in the range m x < 
10 MeV could also generate the correct relic abundance. 
As we have noted, m x ^,5 — 10 MeV is excluded by CMB 
measurements in the standard cosmological scenario (30I ] , 
but smaller masses are possible if the cosmology is mod- 
ified m. 
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FIG. 4: The 
m x , needed 
dark matter 



anapole moment, g/A 2 , as a function of the mass, 
to account for the observed relic abundance of 



IV. DIRECT DETECTION LIMITS 



In this section, we consider the direct detection limits 
on the DM-nucleus scattering cross-section. The Feyn- 
man diagram for this scattering process is shown in Fig. 
[5] According to the derivation in Appendix B, this scat- 
tering cross-section is found to be 
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where mjv and Z are the nuclear mass and charge, v is 
the velocity of x m the lab frame (i.e., the nuclear rest 
frame) and M x n = m x m^ / {m x + m^) is the reduced 
mass of the ^-nucleus system. A similar expression for cr 
has been previously computed in [l[ , but our result differs 
from that expression by a factor of 1 + 2 A4 2 N /m 
ever, if m x <C m^, then M x jv 



and 2 Mj^/m 



2 N . How- 
1| . On the 



reducing our expression to the result of Ref. 
other hand, when to v ~ m^v or m x 3> ton, the fac 



tor 2M 2 N /m 2 N could become significant. In limit where 



i x < ~ '"'iv ^ L X 

m x 3> rriM: we have 1 + 2M xN /m 2 N — > 3. 

A more useful quantity for the direct detection exper- 
iments is 
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where En is the nuclear recoil energy, with a typical value 
~ 1 — 100 keV, and F c (Er) is the nuclear form factor, 
which accounts for the loss of coherent scattering at large 
momentum transfer. 

The differential scattering rate is then 
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where Nt is the number of targets in the detector, p x is 
the local dark matter density, and f(v) gives the velocity 
distribution of the dark matter in the lab frame. The 
lower limit of integration is the minimum velocity needed 
to produce a recoil energy of Er: 
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m N 



y/E R /2m N , 



(11) 



while v m ax is the maximum dark matter velocity in the 
lab frame, given by the sum of the halo escape velocity 
and the velocity of the earth through the halo. 

Eqs. ((9|) and (fT0|) must be convolved with the detec- 
tor efficiency for any given experiment to determine the 
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overall detection rate for dark matter. Here, however, 
we will take a simpler approach (loosely modeled on that 
of Ref. [13), comparing the predicted detection rate for 
anapole dark matter to the detection rate for "standard" 
WIMP dark matter with a nuclear contact interaction 
(CI) . This provides additional insight into the behavior 
of our model compared to WIMP dark matter, but with 
the drawback that the resulting limits on anapole dark 
matter are less precisely determined. 
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More specifically, we calculate the ratio 

(dR/dE R ) a , na . po i e 
(dR/ dEji) c\ ' [iZ) 

where (dR/ dER) &n&vo \ is given by Eqs. ©-([10]), while 
(dR/dEn)ci is the differential scattering rate for a spin- 
independent contact interaction, derived from 



where I\ and Iq are integrals over the dark matter veloc- 
ity distribution given by: 

h = dvf(v)v, (15) 

I = dvf(v)-. (16) 

To compare the limits on WIMP dark matter to the 
limits on anapole dark matter, we select a particular 
value for m x and use the value for g / A 2 shown in Fig. 4. 
For <7 n , we take the boundary of the excluded region for 
a given experiment, and examine the value of TZ as En 
ranges over the values of the recoil energy probed by that 
experiment. Then, if 72. < 1 over this entire range in En, 
the given mass m x with anapole interaction sufficient to 
account for the dark matter is not ruled out by the exper- 
iment under consideration. On the other hand, if TZ > 1 
over the entire Er range, the experimental results do rule 
out this value for m x . In the intermediate regime, where 
TZ changes from < 1 to > 1 over the detectable range in 
En, no conclusion can be drawn without a more detailed 
calculation. 

Eq. (f!4| shows the (Z/A) 2 dependence typical of com- 
parisons between electromagnetically-interacting dark 
matter and dark matter coupling identically to all nu- 
cleons. While Z/A is larger for lighter target nuclei, 
the current XENON100 225 live day results |34| are so 
much more restrictive than any other published limits 
over much of the mass range of interest that we will use 
them here as our standard of comparison. Taking Er 
to lie in the range from 6.6 keV to 43.3 keV, and using 
the truncated Maxwcll-Boltzmann velocity distribution 
for f(v) in Ref. [l2| (ignoring the small seasonal con- 
tribution from the earth's motion around the sun) we 
find that TZ < 1 for all m x < 80 GeV, indicating that 
this mass range for anapole dark matter cannot be ruled 
out by the current XENON100 results. (We have ex- 
amined different models for f(v) and included seasonal 



where A is the total number of nucleons in the target 
nucleus, o n is the dark matter-nucleon cross section, and 
M xn is the reduced mass of the dark matter-nucleon sys- 
tem. Eq. (|13p assumes that the dark matter particle 
couples identically to neutrons and protons. Generally 
speaking, Eqs. (JTUJ) and (TT31 go into the standard cal- 
culation for spin- independent dark matter interactions, 
with the experimental results then used to place an up- 
per bound on o n . 

Using Eqs. ([5]), ([TO ]) , and (fT3" ]l . we find that 



-+(— - J ^)er], (14) 
h \m N 2M 2 N J R y K ' 
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effects to verify that this main conclusion is independent 
of our choice of f(v)). We find that TZ reaches a maxi- 
mum value ^1/4 for m x ^ 30 — 40 GeV, indicating that 
such a particle would lie just below the current thresh- 
old for detection by XENON100. These anapole results 
contrast sharply with the results for electric or magnetic 
dipole dark matter, which allow the lower mass range 
to £ 10 GeV [1,0, HU, but exclude all masses above 10 
GeV. 

Although we restricted our freeze-out calculations to 
m x < 80 GeV, the additional annihilation channels that 
open up at m x ^ 80 GeV will increase the annihila- 
tion rate at a given value of the anapole moment, so 
that the correct relic abundance will be achieved for a 
smaller value of g/ A 2 , with no effect on the DM- nucleus 
cross section. So we expect masses in this higher mass 
range to be even less detectable, and we are justified in 
concluding that m x ^ 80 GeV is also currently allowed 
by XENON100. The main constraint on these heav- 
ier masses comes from the theoretical requirement that 
m x < A, which is violated for sufficiently large masses, 
as is apparent from the calculation in Sec. III. 

V. DISCUSSIONS AND CONCLUSIONS 

Electromagnetically interacting particles provide a rel- 
atively simple model for dark matter: once the magni- 
tude of the interaction is fixed to provide the correct relic 
abundance for a given mass, there is single unique pre- 
diction for the signal in direct detection experiments at 
that mass. (We note in passing that these models can be 
considered an example of isospin-violating dark matter 
proposed by [35[). Unfortunately, for the case of electric 
or magnetic dipole dark matter, the interaction is suffi- 
ciently strong that direct detection experiments already 
rule out such models for the entire range over which such 
experiments have reasonable sensitivity (m ^,10 GeV). 

In contrast, we have shown that a particle interacting 
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exclusively through an anapole moment cannot currently 
be excluded at any mass by direct detection experiments, 
while such a particle would be close to the XENON100 
threshold of detection if m x ^30 — 40 GeV. (Although 
m x ;$ 5 — 10 MeV is excluded by CMB measurements). 
Further, the anapole is the only allowed electromagnetic 
moment for Majorana dark matter. While we have not 
examined collider signatures for this model, the electro- 
magnetic anapole is clearly an interesting new model for 
dark matter and is worthy of more detailed study. 



where a — e 2 /4ir w 1/137 is the fine structure constant 
and we have made the approximation irif -C E « m x . 

The thermally averaged relative velocity in the CM 
frame is given by ^mx)(" r 2 e i) = § T where T is the 
temperature. This implies that (v 2 cl ) = 6T/m x and 
hence 

4q 2 aml ( T \ 

k^//^>=-v^UJ- (A6) 
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Appendix A: The Cross-Section for XX~^ f f 
The scattering amplitude for the process x(p) x(p') ~ * 



M = 



A 2 



u s {k) 7ai tv 0') v r > (p') Y 7 5 u r (p) , (Al) 



where the subscripts r,r' , s, s' are the spin indices of the 
corresponding fcrmions. Averaging over the initial po- 
larizations and summing over the final polarizations, we 
get 



\M\* 



g 2 e 2 
4 A 4 



Tr {(Jf + m/^CJf -m/)7„} 



Tr {(/-m x ) 7 ^ 7 5 (^ + m x ) 7 l ' 7 5 }(.A2) 
In the CM frame, we have p = (E,p), p' = (E,—p), 



k = (E,k) and k' = (E,-k), where E = ^\p\ 2 + m 2 . 
After some algebra, we obtain 



Appendix B: DM-Nucleus Scattering Cross-Section 

The scattering amplitude for the process \(p) N(p') — > 
\(k) N(k') is given by 

M = S «' > MP') M*0 r 7 5 «r(p) , (Bl) 

where the subscripts r,r' ,s,s' are the spin indices of the 
corresponding fermions and Z is the atomic number of 
the nucleus. Averaging over the initial polarizations and 
summing over the final polarizations, we get 

W? = ^^Tr{(^ + m jv )7^(/ + m A r)7,} 

Tr{(^ + m x ) 7 ^7 5 (# + m x )7^7 5 } , (B2) 

where tun is the mass of the nucleus. 

In the CM frame, we have p = (E x ,p), p' = (En, —p), 

k = (E' x ,k) and k' — (E' N ,—k). After some algebra, we 
obtain 
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\M\ 2 = * A4 v, 
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where w rc i = 2v with v = \p\/E being the velocity of 
each annihilating dark matter in the CM frame, and 8 is 
the angle between p and k. 

The differential cross-section is given by 



da xx^ff 



i _ 

1 E 2 



where M x n — m x mN / {m x + mjv) is the reduced mass 
of the x~ nucleus system. Note that \p\ — \k\ — M x n v 
where is v is the velocity of \ in the lab frame. 

In the non-relativistic limit, we have the differential 
cross-section 
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A more relevant quantity for direct detection is 
da/dEn — (dQ/dEn)(da/dn) where Er is the nuclear 
recoil energy. For small momentum transfers, we have 



64 tt 2 (2E) 2 



\M\ 2 . (A4) dn/dE R = 2TTm N /{M 2 N v 2 ) and so 



It is then straightforward to calculate the total cross- 
section 
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Finally, we have the total cross-section 



[1] M. Pospelov and T. ter Veldhuis, Phys. Lett. B 480, 181 
(2000). 

[2] K. Sigurdson, M. Doran, A. Kurylov, R. R. Caldwell, 
and M. Kamionkowski, Phys. Rev. D 70, 083501 (2004); 
erratum, Phys. Rev. D 73, 089903 (2006). 

[3] S. Gardner, Phys. Rev. D 79, 055007 (2009). 

[4] E. Masso, S. Mohanty, and S. Rao, Phys. Rev. D 80, 
036009 (2009). 

[5] A. L. Fitzpatrick and K. M. Zurek, Phys. Rev. D 82, 
075004 (2010). 

[6] W. S. Cho, et al., Phys. Lett. B 687, 6 (2010); erratum, 

Phys. Lett. B 694, 496 (2011). 
[7] J. H. Heo, Phys. Lett. B 693, 255 (2010). 

[8] J. H. Heo, Phys. Lett. B 702, 205 (2011). 

[9] T. Banks, J.-F. Fortin, and S. Thomas. larXiv:1007.55l"5l 

[hep-ph]. 

[10] V. Barger, W.-Y. Keung, and D. Marfatia, Phys. Lett. B 
696, 74 (2011). 

[11] J.-F. Fortin and T. M. P. Tait, Phys. Rev. D 85, 063506 

(2012). 

[12] E. Del Nobile, et al.. larXiv: 1203.6652 1 [hep-ph], 

[13] V. Barger, W.-Y . Keung, D. Marfatia, and P.-Y. Tseng, 

larXiv:1206.0640l [hep-ph] . 

[14] J. H. Heo and C. S. Kim. larXiv: 1207. 1341 II [astro-ph] . 
[15] Ya. B. Zel'dovich, Sov. Phys. JETP 6, 1184 (1958). 
[16] V. M. Dubovik and V. E. Kuznetsov, Int. J. Mod. Phys. 

A 13, 5257 (1998). 
[17] C. S. Wood, et al., Science 275, 1759 (1997). 



[18] A. L. Fitzpatrick and K. M. Zurek, Phys. Rev. D 82, 

075004 (2010). 
[19] E. E. Radescu, Phys. Rev. D 32, 1266 (1985). 
[20] F. Boudjema, C. Hamzaoui, V. Rahal and H. C. Ren, 

Phys. Rev. Lett. 62, 852 (1989). 
[21] E. N. Bukina, V. M. Dubovik and V. E. Kuznetsov, Phys. 

Lett. B 435, 134 (1998). 
[22] B. Holdom, Phys. Lett. B 166, 196 (1986). 
[23] N. Arkani-Hamed, D. P. Finkbeiner, T. R. Slatyer and 

N. Weiner, Phys. Rev. D 79, 015014 (2009). 
[24] J. D. Bjorken, R. Essig, P. Schuster and N. Toro, Phys. 

Rev. D 80, 075018 (2009). 
[25] M. Pospelov, Phys. Rev. D 80, 095002 (2009). 
[26] E. W. Kolb and M. S. Turner, The Early Universe, (New 

York: Addison- Wesley, 1990). 
[27] R. J. Scherrer and M. S. Turner, Phys. Rev. D 33, 1585 

(1986). 

[28] G. Steigman, B. Dasgupta, and J. F. Beacom, Phys. Rev. 

D 86, 023506 (2012). 

[29] C. M. Bender and S. Sarkar. larXiv:1203.1822l [hep-th]. 

[30] C. M. Ho and R. J. Scherrer. larXiv: 1208.43471 [astro-ph], 

[31] K. Griest and D. Seckel, Phys. Rev. D 43, 3191 (1991). 

[32] E. Komatsu, et al., Astrophys. J. Suppl. 192, 18 (2011). 

[33] C. M. Ho and R. J. Scherrer, in preparation. 

[34] E. Aprile, et al.. larXiv: 1207.5988 1 [astro-ph], 

[35] J. L. Feng, J. Kumar, D. Marfatia and D. Sanford, Phys. 

Lett. B 703, 124 (2011). 



